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Existence of strong solutions in a larger space for the 

shallow-water system 

Frederic Charve* Boris Haspot ^ 



Abstract 



This paper is dedicated to the study of both viscous compressible barotropic fluids 
fvj ' and Navier-Stokes equation with dependent density, when the viscosity coefficients 

are variable, in dimension d > 2. We aim at proving the local and global well- 
posedness for respectively large and small initial data having critical Besov regularity 
and more precisely we are interested in extending the class of initial data velocity 
when we consider the shallow water system, improving the results in [fl, 24] and 
[22]. Our result relies on the fact that the velocity u can be written as the sum 
i-G ' of the solution ul of the associated linear system and a remainder velocity term 

C^ , u; then in the specific case of the shallow-water system the remainder term u is 

more regular than ul by taking into account the regularizing effects induced on the 
bilinear convection term. In particular we are able to deal with initial velocity in 
iJ2^i as Fujita and Kato for the incompressible Navier-Stokes equations (see [19]) 
K^ , with an additional condition of type mo G B^ i. We would like to point out that 

fvj ' this type of result is of particular interest when we want to deal with the problem of 

^T) , the convergence of the solution of compressible system to the incompressible system 

CO ■ when the Mach number goes to 0. 

(N 

: 1 Introduction 

The motion of a general barotropic compressible fluid is described by the following system: 

' dtp + div{pu) = 0, 

< dt{pu)+dw{pu'^u)-div{n{p)D{u)) -V{\{p)divu) + VP{p) = pf, (1.1) 

^ {p,u)/t=o = {po,uo). 

Here u = u(t, x) G W^ stands for the velocity field and p = p{t, x) G M+ is the density. 
The pressure P is a suitable smooth function depending on the density p. We denote 
by A and p the two viscosity coefficients of the fluid, which are also assumed to depend 
on the density and which verify some parabolic conditions for the momentum equation 
p > and X + 2p > [in the physical cases the viscosity coefficients verify X + -^ > 
which is a particular case of the previous assumption) . We supplement the problem with 
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initial condition {po,uq) and an external force /. Throughout the paper, we assume that 
the space variable x G W^ or to the periodic box T'^ with period Oj, in the i-th direction. 
We restrict ourselves to the case d >2. Let us recall that in the case of constant viscosity 
coefficients, existence and uniqueness for (1.1) in the case of smooth data with no vacuum 
has been stated in the pioneering works by Nash (see [29]) and Matsumura and Nishida 
( see [25, 26]). 

In this article we obtain the existence of strong solution (in finite time with large initial 
data or in global time for small initial data) for an optimal class of initial data by combin- 
ing two different ingredients, first the notion of scaling and seconds by taking advantage 
of suitable choices on the viscosity coefficients which may confer specific structures in 
terms of regularity. We will detail more this last point later, but as a first example of the 
importance of the viscosity coefficients we recall that in [27] A. Mellet and A. Vasseur 
have obtained the stability of the global weak solutions for the Saint- Venant system by 
using new entropy giving additional regularity on the gradient of the density and on the 
integrability of the velocity. 

Let us recall the fundamental notion of scaling for system (1.1). Indeed guided in our 
approach by numerous works dedicated to the incompressible Navier-Stokes equation (see 
in particular the pioneering works of Fujita and Kato concerning the existence of strong 
solutions for the incompressible Navier-Stokes equations in [19]), we aim at solving (1.1) 
in the case where the data {pQ,uo) have critical regularity. By critical, we mean that we 
want to solve the system (1.1) in functional spaces with invariant norm by the natural 
changes of scales which leave (1.1) invariant. More precisely in the case of barotropic 
fiuids, the following transformations: 

{p{t, x),u{t, x)) — > {p{l^t, Ix), lu{l^t, Ix)), I G M, (1.2) 

verify this property, provided that the pressure term has been changed accordingly. This 
notion of critical functional frameworks has been extensively used in order to obtain opti- 
mal class of initial data for the existence of global strong solution (see [6, 10, 13, 18, 21, 24] 
in the case of constant viscosity coefficients). In particular the first result on the exis- 
tence of strong solutions in spaces invariant for the scaling of the equations (when the 
viscosity coefficients are constant) is due to R. Danchin in [13] when the initial data 

_ _d _ _d -I 

{qo = Po ~ P^fJ-o) (with p > 0) are in B21 x {B21 )"'. In [18], R. Danchin generalizes the 

d d 1 

previous results by working with more general Besov space of the type B^^^ x {BI!^^ ) 
with some restrictions on the choice oip {p < d) due to some limitation in the application 
of the paraproduct law. The fact that p = pi \s a. consequence of the strong coupling 
between the density and the velocity equations, indeed the pressure term is considered 
as a remainder for the parabolic operator in the momentum equation of (1.2). In [24], 
the second author generalizes the results of [18] inasmuch as we have no restriction on 
the size of p for the initial density. To do this we are working with a new variable, the 
effective velocity which allows to cancel out the coupling between the pressure and the 
velocity. 

In the case of global strong solution in critical space for small initial data, we would 
like to recall the works of R. Danchin in [14] who shows for the first time a result of 
global existence of strong solution close to a stable equilibrium when the initial data 

verify (qoj^o) £ (-^2^1 '^ -^2^1 ) ^ -^21 • "^^^ main difficulty is to get estimates on the 



linearized system given that the velocity and the density are coupled via the pressure. 
What is crucial in this work is the smoothing effect on the velocity and a L^ decay on 
p — p (this plays a key role to control the pressure term). This work was generalized in 
the framework of Besov space with a Lebesgue index different of p = 2 by Q. Chen et al 
in [10] , F. Charve and R. Danchin in [6] and B. Haspot in [22]. 

However very few articles really take into account the structure of the viscosity coef- 
ficients, indeed most of them deal with constant viscosity coefficients or consider the 
system (1.1) as a perturbation of the previous case. In addition to have a norm invari- 
ant by (1.2), appropriate functional spaces for solving (1.1) must provide a control on 
the L°° norm of the density for at least two reasons; the first corresponds to avoid the 
vacuum and to ensure the parabolicity of system (1.1), the second is to assume some 
property of multiplier on the density in order to deal with the term -A. It explains 
why the authors restrict their study to the case where the initial density is assumed in 

d d 

Bpi with p G [l,+oo[ suitably chosen, indeed Bp^ is embedded in L°° . Furthermore 
in order to propagate this regularity on the density via the mass equation, it appears 

d 

necessary to assume that the velocity is Lipschitz, it means Vn G L\,{B^i) ^-J- L}p{L'^) 
with p G [l,+oo[ and T > 0. We would like to point out that this necessary Lipschitz 

control on the velocity seems to prevent any hope of existence of strong solutions when 

■ - — 1 
uq is only assumed to belong in H2 

Indeed in [6, 10, 13, 18, 21, 24] the idea is to propagate a L°°{Bl'^^) n L^Bl^^) regu- 
larity on the velocity u via the regularizing effects induced by the momentum equation 
written in its eulerian form: 

dtu - -An - ^^^^Vdivu = R, (1.3) 

P P 

with: R = —u • \7u -^ + /. However if we only were interested in obtaining a 

regularity in L°°{H^^ ) n L^{H^^ ) , it would be necessary to deal with the term -An 
and getting enough regularity in order that - remains in a multiplier space M{H^^ ) 

of H^^ . Typically H^ Pi L°^ is embedded in Ai{H^^ ). However in this case how 

■ i 
to propagate the regularity H2 f] L°° on the density as our velocity is only assumed in 

L^{H^^ ) (it means not necessary Lipshitz)? 

We want to partially solve this question in the case of specific viscosity coefficients. As 

we explained above, one of the main difficulty is linked to the treatment of heat equation 

with variable coefficients. We would like to work with the shallow water system (i.e 

p{p) = HP and A(p) = 0, and to simplify we will take p = 1) that we can rewrite in the 

following form: 

' dtp + div(pu) = 0, 

< dtu + u-Vu- Au-2D{u).V {In p) + V{G{p))=0, (1.4) 

, {p,u)/t=o = {po,uo), 



where the operator A is defined by: Au = Au + Vdivu and G (p) = — —. Roughly 



Pjp) 
p 
speaking, we are led to consider a basic heat equation with remainder terms, and in par- 
ticular it seems possible to propagate in this case the regularity L°^{H2~ ) n L^{H2^ ). 



However due to the strict coupling between the pressure and the velocity, it will be nec- 

d 

essary to control the density in L}p{B 21) (the third index 1 is due to multiplier space 

reason). More precisely we will split the solution u into the following sum u = ul + u 

with: 

f dtUL - AuL = 

\ (1-5) 

[ w/t=o = ^0 

In the present article, if we assume that uq belongs to H^~ n B^ ^ then, with usual 

■ — +1 
methods we show that u is more regular than ul {u will be in L}j.{B2i )) which will help 

us propagating the regularity on u. 

To simplify the notation, we assume from now on that p = 1. Hence as long as p does 

not vanish, the equations for [q = p — \^u) read: 

{dtq + u-Vq + {l + g)divu = 0, 

\dtu + u-Vu-Au- 2D(u).V(ln(l + q)) + V(G(1 + g)) = 0, 

We can now state our main result: 

Theorem 1 Let P he a suitably smooth function of the density. Assume that uq G 

d_ 1 _ _ d o , d -1 

-62^2 ^ ^00 1' divuo G B21 , qo € B21 and that there exists c > such that po > c. 
Then there exists a time T such that there exists a unique solution {q, u) for system (1-1) 
on [0, T] with 1 + q bounded away from zero and, 

qGCriBl,) and u G CTiBl2') ^i L^risl^^) nCriB^],) n L'^bI^i). 

— - — — I ' —4-1 

In addition divu belongs to L^{B2i ) Pi L^(i?2^^ ). 

Moreover if P'{1) > 0, there exists a constant eq such that if: 

Ikoll . d_i d + \\uo\\ . d_i < £0, 

then the solution is global. 

Remark 1 We would like to mention that this theorem could easily be extend to the 
case of Besov space constructed on Lebesgue spaces with index p / 2. More precisely as 
in [13], we could obtain the existence of solution under the condition than 1 < p < 2N 

d_ d 1 

and the uniqueness for I < p < N when {qo,uo) belong to B^^ x {Bp2 H B^ ^). These 
restrictions on the size of p are essentially due to some limitations on the use of the 
paraproduct when we are dealing with the convection term u ■ Vu. 

Remark 2 In [5] the authors obtain a new entropy inequality for the shallow-water 
system. They show that we can control ^ in L°°[L'^). Roughly speaking it means that 
we control the density p in L°°{H^). It means in particular that our initial data are 
very close to the energy data in dimension d = 2. Indeed we essentially assume only 
additional condition in terms of vacuum ( — G L°°) and a control on the L°° norm of 
Pq. This conditions are natural in order to deal with the non linear terms but also for 



preserving the parabolicity of the momentum equation. Moreover we suppose that uq 
is in B^ ^. We would hke to emphasize that this last condition is quite optimal for 
Navier-Stokes system. Indeed in the case of the incompressible Navier-Stokes equations 
Bourgain and Pavlovic have proved that the system is ill-posed when uq £ B^oo f^ee 

[4])- 

Remark 3 This result is also very interesting in the case of the convergence of the 

solution of compressible system to the incompressible system when the Mach number 

goes to 0. Indeed we are able to deal with initial data of Fujita-Kato type, it means 

here uq G i^2^ n B^ ^ which improves the result of R.Danchin (see [16]) who needs to 

■ i-l 
assume uq G B^ i . 

d rx 

Remark 4 The additionnal assumption divuo € B21 is quite natural as it expresses 
the compressibility of the Euid. 

The previous theorem can be easily adapted to the incompressible density dependent 
Navier-Stokes equations. We recall here the equations: 

' dtp + div{pu) = 0, 

dt{pu) + diY{pu(g)u) -div{2pDu) + \7U = pf, 

divu = 0, 
, {p,u)/t=o = {po,uo). 
Following the same ideas than in theorem 1, we obtain the following result: 

Theorem 2 Let d > 3. Assume that uq G -62^2 ^ ^00 1 '^^^h divuo = 0, qq £ i?! 1 and 
that there exists c > such that po > c. Then there exists a time T such that there 
exists a unique solution {q,u) for system (1.1) on [0,T] with 1 + q bounded away from 
zero and, 

qeCriBh) and uGCTiBl,')nL},{Bl^')nCTiB^\)nL},iBl,) 



and Vne L^{B^2^] 



Moreover there exists a constant Eq such that if: 



|<?o|| . 4 + \Wo\\ . d_^ < Eq, 



^2,1 -°2,2 



then the solution is global. 



Remark 5 This result may be considered as an extension of the Fujita-Kato theorem 
(see [19]) to the case of incompressible density dependent Navier-Stokes equations. In 
particular it improves the analysis of [1] and [23] inasmuch as we deal with a velocity in 
a Besov space such that the third index is different of 1 and with a critical initial density 
in terms of scaling. 

Our paper is structured as follows. In section 2, we give a few notation and briefly 
introduce the basic Fourier analysis techniques needed to prove our result. In section 
3 and section 4, we prove theorem 1 and more particular the existence of such solution 
in section 3 and the uniqueness in section 4. In section 5 we are proving the global 
well-posedness of theorem 1. In section 6, we are dealing with theorem 2. 



2 Littlewood-Paley theory and Besov spaces 

As usual, the Fourier transform of u with respect to the space variable will be denoted 
by J-{u) or u. In this section we will state classical definitions and properties concerning 
the homogeneous dyadic decomposition with respect to the Fourier variable. We will 
recall some classical results and we refer to [2] (Chapter 2) for proofs (and more general 
properties). 

To build the Littlewood-Paley decomposition, we need to fix a smooth radial function 
X supported in (for example) the ball i?(0, 3), equal to 1 in a neighborhood of i?(0, |) and 
such that r 1— )■ xif-^r) is nonincreasing over IR-|_. So that if we define (p{$,) = xi^/'^)~xiO^ 
then if is compactly supported in the annulus {.^ G M , | < |^| < 1} and we have that, 

V^GM'^MO}, J]99(2-'0 = l. (2.8) 

lez 

Then we can define the dyadic blocks {iS.i)i^i by A; := Lp{2^^D) (that is A^n = ip{2~^ ^)u{$^)) 
so that, formally, we have 

n = ^AiU (2.9) 

I 

As (2.8) is satisfied for ^ 7^ 0, the previous formal equality holds true for tempered 
distributions modulo polynomials. A way to avoid working modulo polynomials is to 
consider the set S'^ of tempered distributions u such that 

lim llS'^nllioo = 0, 
i— >— 00 

where Si stands for the low frequency cut-off defined by Si := x(2~'-C'). If n G 5^, (2.9) 
is true and we can write that Siu = > A^n. We can now define the homogeneous 

k<l-l 

Besov spaces used in this article: 

Definition 1 For s G M and 1 < p,r < 00, we set 



\u\\b3 



y^2'"'''||Ain||£p ) if r < 00 and \\u\\r,s := sup2''*||A/u|kp. 



We then define the space B:Lj. as the subset of distributions u ^ S'^ such that ||n||^s is 
finite. 

Once more, we refer to [2] (chapter 2) for properties of the inhomogeneous and homoge- 
neous Besov spaces. Among these properties, let us mention: 

• for any p G [1, 00] we have the following chain of continuous embeddings: 

Bl, ^ LP ^ Bl^- 

d 

• if p < 00 then BI^ ^ is an algebra continuously embedded in the set of continuous 
functions decaying to at infinity; 



for any smooth homogeneous of degree m function F on M \{0} the operator F{D) 
maps ijp,, in B^""^. This imphes that the gradient operator maps B^^^ in B^^^. 



We refer to [2] (lemma 2.1) for the Bernstein lemma (describing how derivatives act on 
spectrally localized functions), that entails the following embedding result: 

Proposition 1 For all s £ M, 1 < pi < p2 < oo and 1 < ri < r2 < oo, the space -B^^ ^.^ 



' embedded m the space "" 
Then we have: 



. s — ^y ^ 

is continuously embedded in the space Bp2,r2^^ ^^ 



In this paper, we shall mainly work with functions or distributions depending on both the 
time variable t and the space variable x. We shall denote by C{I] X) the set of continuous 
functions on / with values in X. For p G [1, oo], the notation ^^(1; X) stands for the set 
of measurable functions on / with values in X such that 1 1— t- ||/(t)||x belongs to LP(I). 

In the case where I = [0,T], the space LP{[0,T];X) (resp. c([0,T];X)) wih also be 
denoted by Lj,X (resp. CtX). Finally, if / = M+ we shall alternately use the notation 
LPX. 

The Littlewood-Paley decomposition enables us to work with spectrally localized 
(hence smooth) functions rather than with rough objects. We naturally obtain bounds for 
each dyadic block in spaces of type Lf^LP. Going from those type of bounds to estimates 
in Lj^Bpj. requires to perform a summation in €"{1,). When doing so however, we do 
not bound the Lj,Bp ^ norm for the time integration has been performed before the F' 
summation. This leads to the following notation (after J.-Y. Chemin and N. Lerner in 
[9]): 

Definition 2 For T > 0, s G M and 1 < r, p < cx), we set 

II^IIl^b|,, •- Ir ll^<?^llL^Lp|lr(z)- 

One can then define the space Lf^B^j. as the set of tempered distributions u over (0, T) x 
M'^' such that \imq^_ooS„u = in LP([0,T];L°°(R'^)) and IkllTp a, < oo. The letter T 

is omitted for functions defined over IR+. The spaces Lf^B^^ may be compared with the 
spaces L!^Bp ^ through the Minkowski inequality: we have 

IkllrpR^ <||'u||rPRs if r > /9 and ||n||TpA^ >\\u\\rPns if r < p. 

^T P,r -^T P,r ^T P,r ^T P,r 

All the properties of continuity for the product and composition which are true in Besov 
spaces remain true in the above spaces. The time exponent just behaves according to 
Holder's inequality. 

Let us now recall a few nonlinear estimates in Besov spaces. Formally, any product 
of two distributions u and v may be decomposed into 

uv = TuV + T^u + R{u,v), where (2.10) 

TuV := 2^ Si-iuKiv , T^u :=y^ Si-ivAiu and R{u,v):=y^ \^ AiuAi'V. 
I I I |/'-«|<i 



The above operator T is called "paraproduct" whereas R is called "remainder". The 
decomposition (2.10) has been introduced by J.-M. Bony in [3]. 

In this article we will frequently use the following estimates (we refer to [2] section 
2.6, [14], for general statements, more properties of continuity for the paraproduct and 
remainder operators, sometimes adapted to Lj,Bp ^ spaces): under the same assumptions 
there exists a constant C > such that if l/pi + l/p2 = ^/p, and 1/ri + l/r2 = 1/r: 

IIT'u'vIIbs < C'll'"llL°°ll'y||Rs , 

\\Riu,v)\\.s,+s,-i<C\\u\\^s, \\v\\j^s, (si + s2>0). (2.11) 

Let us now turn to the composition estimates. We refer for example to [2] (Theorem 
2.59, corollary 2.63)): 

Proposition 2 l.Lets>0,ue B^^ D L°° and F G w/„'J+^'°°(M'^) such that F(0) = 
0. Then F{u) G -Bf i ^^^ there exists a function of one variable Cq only depending 
on s, d and F such that 

\\F{u)\\^s^^<Co{\\u\\l^)\\u\\^^^^. 

2. Ifu and v G S|i and if v - u £ B^^ for s g] - f , f] and G £ wl^l^^'°°{R'^), then 
G{v) — G{u) belongs to i?| ^ and there exists a function of two variables C only 
depending on s, d and G such that 

\\G{v) - G{u)\\^.^ ^ < C{\\u\\l^M\l^) { \G'm + \\u\\ . . + \\v\\ .A\\v- ^b| ,• 

\ ^2,1 ^2,1/ 

Let us now recall a result of interpolation which explains the link between the space Bt^ ^ 
and the space Bt,^ (see [12] or [2] sections 2.11 and 10.2.4): 

Proposition 3 There exists a constant C such that for aii s G M, e > and 1 < p < +oo, 

II II ^ /^ T^ II II 111 ^T^ P^°°) "l^j'KOpaa) 

\\u\\iP,^s -, < C; \\u\\iP,^s ^log e + 

^Ty^p,oo) 

Let us end this section by recalling the following estimates for the heat equation: 

Proposition 4 Let s G M, (p, r) G [l,+oo]^ and 1 < /C2 < Pi < +oo. Assume that 
uq G Bpj. and f G L^{Bp^r '^). Let u he a solution of: 

J dtu - iJ.Au = f 

\ Ut=0 = UQ . 

Then there exists G > depending only on N, fi, pi and p2 such that: 
If in addition r is finite then u belongs to C([0, T], B^^^). 



3 Existence of solution 

3.1 A priori estimates 

Let us emphasize that we do not use the fact that P'{1) > 0. Denoting by G the unique 
primitive of x i— t- P'{x)/x such that G{1) = 0, recall that system (1.6) now reads: 

{dtq + u-Vq+{l + g)divu = 0, 

\dtu + u-Vu-Au- 2D(u).V(ln(l + q)) + V(G(1 + ^)) = 0, 

where the operator A is defined by: Au = Au + Vdiv u. 

Let ul be the unique global solution of the following linear heat equation: 

(d,u,-Au, = 0, ^^^^^ 

[ ULt=0 = ""0, 

Thanks to the classical heat estimates recalled in (4) (we refer for example to [2], lemma 

■ --1 
2.4 and chapter 3), as no G -62^2 ^ ^00 1 "^^ have for all time: 

UL G (irBl-'nLlBl;^') n (LrB^\nLlBl^,y (3.14) 

and the corresponding energy estimates. Moreover, as divuo G B21 , and as divu^ 

satisfies: 

f Stdiv UL - AuL = dtdiY ul - 2Aul = 0, /„ -, ,-^ 

< (3.15) 

divuL|t=o = divtio, 



we also have that: 



divuL G L^B^^^ n Lj^li, (3.16) 



which will be crucial in the study of the density equation. Then, if we denote by u = 
u — Ul, we now need to study the following system: 

' dtq +(u + UL).Vq + (1 + ^)div {u + ul) = 0, 
dtu — Au — Vdiv n + (n + ul)-Vu + u.Vul + ul-Vul 

- 2D{u + nL).v(ln(l + g)) + v(g(1 + ^)) = 0, 

(3.17) 
The interest of introducing this system is that the most problematic term in the addition- 
nal external force terms, namely ul ■ ^ul, is in fact regular. Thanks to the paraproduct 
and remainder estimates, we have: 

\\ul-^ul\\ ^ .4_i < IKlIIjoc^-i IIulII d+i 

In this article we will prove existence and uniqueness of a local solution such that the 

~ ■ ——1 • — +1 

velocity fluctuation u is in the space L^i?! 1 H L}pB2i for some T > 0. Then, going 

back to the original functions, thanks to the following embeddings: 

B~C-^B-l2\ bI-'^B^„ 



we will end with a velocity u with the same regularity as ul (that is 3.14 and 3.16). 

Let us now state the following transport-diffusion estimates which are adaptations of 
the ones given in [2] section 3: 

— -1-1 

Lemma 1 Let T > 0, -| < s < |, uq G -Bf,!, / G -^t-^I.i ^^'^ v,w e -^^-^1,2 ^ ^L,i- 
If u is a solution of: 

{dfU + V ■ Vn + u ■ Vw — Au = f, 
Ut=o = Uo, 

then, ifV(t) = L (||Vt;(r)|| d + 11 Vttjfr) 11 d jdr, there exists a constant 



C > such that for ah t G [0, T]; 



Proof: We refer to [2] for details. Localizing in frequency, if for j G Z, Uj = AjU, we 
have: 

dfUj + V • Vuj — Auj = fj — Aj{u • Vw) + Rj, 

where Rj = [v.\7, Aj]u. Taking the L^ innerproduct we obtain: 

9t\\uj\\L2 +2^-'||tij||L2 < ||divi;||Loo||nj||j^2 + ||/j||l2 + \\Aj{u- Vw)\\l2 + ||it'j||L2. 

Classical commutator estimates (we refer to [2] section 2.10) then imply that there exists 
a summable positive sequence Cj = Cj{t) whose sum is 1 such that: 

||i?J| r2 < Co2~-^*||Vf II d ll^tllns . 

JUL - J 'bI^HBI/ "^M 

Thanks to the paraproduct and remainder laws, we have: 

||ti-Vtf||Rs < C||n||rjs ||Vw|| d 

11^2,1 - 11^2,1 II "sl.nBSo,! 

SO that we finally obtain the result. ■ 

Concerning the transport equation of the density fiuctuation, localization implies: 
dtQj + {u + ul) ■ Vqj = -Aj ({1 + g')div {u + ul)) + Rj, 
where Rj = [{u + ul)-V, Aj]q. Using the same method, we get the estimate: 

|k|L d < IkolU d + / (||g|| d ||V(li-|- ul)|| d 

L^Bl, L^Bl, Jo V Bl, Bl,nBO^, 

+ (1 + Ikll .| )(l|divuL|| .d +\\Vu\\ ,d))dT. (3.18) 

Bli B^A B^A ' 
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Remark 6 Note that here, due to the external force terms for the density, we need to 

d 

have a B2 i-control ofdivui- 
With a rough majoration, 






,1 -^0 ■02,1 ^ 



where 



V(t) = / I ll^^ll d,-i + llVurll d + lldivtir II d ]dT, 

/n V B^ B^ nfiO B^ J 

•'U ^2,1 -°2,2' '^00,1 ^2,1 

and then thanks to the Gronwah lemma, 

IklU . <Ce'=^^W(|koL . +/V'(r)e~^^Wdr), 

which gives: 

Ikll d < Ce^^(*M l + |ko|| d 1-1, (3.19) 

For the velocity, the same localization technique is used, and if we introduce: for j € TL, 
Uj = AjU, we have: 

dtUj + {u + ul) • Vuj — Auj = Ajf — Aj{u • Vu^) + Rj, 

where 

(Rj = [{u + ul).V,Aj]u, 

\ f = -UL.VUL + 2D{u + UL).v(ln{l + q)) - v(g(1 + < 



Taking the L^ innerproduct with Uj we obtain: 

9t\\Uj\\L'2+2^^\Uj\\L2 < ||div(lZ+UL)||L^||Mj||^2 + ||Aj/||j^2 + \\Aj(u-Vul)\\l2 + ||-Rj||l2- 

After a time integration, the mutliplication by 2-'^2~ •* followed by a summation over 
J G Z gives {u{0) = and the commutator is estimated as above): 



roor>5 rlo^ 

^t -°2,1 ^t-°2,l 



d , + \\u\\ d , , < 

ft 

'0 -02,1 -°2,1 -°2,2' '-°oo,l -°2,1 

Classical computations show that: 



f (\\f\\.d, + \\u-VuL\\.d, + \\Viu + UL)\\d \\u\\.d,)dT. (3.20) 

In \ R" R" B" nR" R" / 

./U -°2.i -°2.i -"2 -;■' '-°no 1 -°2.1 



llti-Vtiill d_-, < ||u|| d_n||VnLll d , 

and from the definition of /: 

. d, < \\ul\\b-\ \Wl\\ . 4+1 + ||VG(1 + g)|| . d, + 2||Vln(l + q) • D{u + nL)|| . d,. 

^2,1 ^2,2 -"2,1 ^2,1 
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The second term is estimated thanks to the composition lemma (see ???): 



||VG(1 + (?)||.. , <||G(1 + (7)||.. <Ci\\q\\L^: 



B4 



B4 



The last term has to be treated carefully. As we cannot rely on the smallness of q (our 
data can be large), we will follow the same method as R. Danchin in [18] and thanks to 
a frequency cut-off we decompose this term into two parts which will both be small. For 

771 G Z: 



I = V(ln{l + q)-ln{l + Smq)) -Diu + UL) 



Vln{l + q)-D(u + UL) = I + II with: 

[// = V ln(l + Smq) ■ D{u + ul). 

Thanks to the paraproduct and remainder laws we have, 



/|| .d_i < ||ln(l + g)-ln(l + S'„g)|| .d \\V{u + ul) 



^2,1 



Bl2^B'L,i 



<C{\\q\\L^)\\q-S„,q\\, \\V{u + Ul)\\, 



B4 



B„\nB' 



■0 



(3.21) 



We then use the estimate given in [18]: there exists a constant C > 1 such that: 
\\q - Sn,q\\ . d < \\qo - SmqoW -i + (1 + Ikoll . 4 ){e^'^^'^ - 1) 

^2,1 -t>2,l B^-^ 

If < a < 1, the other term is estimated the following way: 



\II\\^4-i < l|Vln(l + 5„g)||_^4+_J|D(tI + nL)||^d_ 



'B^~ 
^2,1 



B^' 

-°2,1 



<Ci\\q\\L^)2"'''\\q\\. ||n + ^L||.4+i-.- (3-22 



-°2.1 



B^ 

-°2.2 



And thanks to the Gronwall lemma, we finally obtain that for all < a < 1 (for example 
we can take a = 1/2) and all ttt, G Z: 






d 1 + \\u\\ 



< 



L]B^J 



^^L TooR-i \\Ul 



t^2,2 •'O y «2,1 ^2,2 



+ IIq'o - SjnqoW d + (1 + Ikoll .d){e 



CV{r) 



-°2,1 



-°2,1 



l))||u + nL| 



-°2,2 ' '-°oo,l 



dT 



(3.23) 



Thanks to the Holder estimate, we have: 



\u + ul\\ , .d_^i_„ < t2 U + Ul _i 

lt^bI:'- 



and by interpolation, 



a 



< 



a 



B4 



LT^l^ L\B? 



'w < - 
--^1 - 2"-' "Joo^f 
2,2 * 2,2 



-x+(l 



a 



-'^t-°2,2 
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so that, 



u + 



r 1 o2 TOO n^ 



If we denote /?(t) = llulL d_, + ||ii|| d, , , the estimates on the velocity turn into: 



/3(t)<e^^W llnillroco-i 



PL I 



IR^-I 



^t'^: 



L°°B^^ 



L?°B4 



+C{\\qL .4) ^IklU .4+ ll'?o-5™go||.4+(l+lko||.4)(e^'^*^-l)+2"'"t-|klL 



Bo^ 



■'"t -"2,1 



X /?(t) + ||tXL|U .d , . ^ +||nL|| ,^, . (3.24) 

Let us now state and prove the fohowing lemma: 

Lemma 2 Let {q,u) satisfying (SW) on [0,T] x R'^. Assume that q G C^{[0,T],Bl^) 
and u, UL G C^{[^, T],B^2 ^ -^oo\ ^ -^22 ^ -^i. i))"*' where ul satisfies: 



dtUL-AuL = 0, UL\t=0 = Uo. 



(3.25) 



If we denote u = u — ul, there exist three positive constants rj, C (only depending on d), 
C (depending on C, d and qq) and m £ Z (= m{ri, qo)) such that if rj g]0, 1] satisfies: 



rje 



3^''(1 + C")(1 + 



\Uo\\ . d_ 



< 



1 



-^2,2 ri-B^^ 1 



(3.26) 



and m is chosen such that \\qo — SmqoW '^ ^ v'^- Then ifT is small enough so that: 

-°2,1 



/n (l|VtiL||.d +||divnL|| .4 )c?r< 77^ 

TC + (1 + C")(e^^(^) - 1) + 2°™TtC" < if. 
then we have, for all t G [0, T], 



(3.27) 



\q\\~ ■'i 

Loo d"2" 
rri -Do 1 



< e^^^' 1 + 



^0 d 

Bh 



L?°B. 



d , + \\u\\ d,, < 2ri. 



(3.28) 



Proof: We recall that {q, u) satisfies on the interval [0, T] the following system: 

' dtq +{u + UL).Vq + (1 + g)div {u + ul) = 0, 
dfU — Au — Vdivn + (n + ul).Vu + u.Vul + ul.Vul 

- 2D{u + nL).v(ln(l + g)) + v(g(1 + ^)) = 0, 
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Assume that T is small enough so that we have for some r/ G [0, 1] (to be precised later): 







II ViiL II. d . +||divuL|| d ]dT <r]^ <ri, 



and let us define 



r* = sup{tG [0,T], j ||n|Ld+i <2r?}. 
Then for ah t G [0,r*[, from (3.19) we have (r/ < 1): 



\q\\_ d 



< ^C(2r,W) 1 + 



^0 d 

b7. 



def 



l<e^^ 1 + lkoll^d -1 = CO. (3.29) 



Concerning the velocity, using the previous estimate on q and the fact that 



kill- i-i . , +II^lII , 



< 



■^(-^2,2 ri-^oo.i 



allows us to write that: 



=,3Cr, 



Poll . d_i 

-^2,2 ri-B^ 1 



CV(t). 



/3(t) <e^'-'' ll^olls-i r + C(co) teo+ Iko-^mgoll .d +(l+co)(e^^^^^-l) + 2"'"t2co 



X (/3(t) + ||no||.d , . Jl ). (3.30) 
V R,?„ nB^^/J / 



s|,2 ni?;,;i^ 

Let us first fix m = m(ry) such that ||go — 5'm<?o|| d < ??^- Then let us take T so small 
that: 

The estimate turns into: 



Tco + (1 + co)(e^^(^) - 1) + 2"'"rtco < ry^. 



/3(t)<e3^M |hollB-i^' + C(co)^'[l + /3(t) + ||no||.d , . ^ 

\ oc,l L R_2_ nR ^ 

Then, if r/ is taken so small that 

^e^''''(ho|lB-i +C(co)(l + ||t/o||.d , . J< 

oo,l R?. nR"^. 

then as 77 g]0, 1], we obtain that: 



Bl.'nB-- 2' 



m<^m)+ri), 



which implies: 



m < V- 

Then by contradiction this leads to T* = T and for all t £ [0, T], 

ll'ull d , + ||u|| d , 1 < 2?7. 

^t ^2,1 ^t^2,l 
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3.2 Existence 

We use a standard scheme for proving the existence of the solutions: 

1. We smooth out the data and get a sequence of smooth solutions (g", u^)nef^ to (1.6) 
on a bounded interval [0, T"] which may depend on n. 

2. We exhibit a positive lower bound T for T", and prove uniform estimates on (g", u") 
(we refer to the next subsection for the definition of u") in the space 

Et = Ct{bI,) X (CT(i?|i)n4(i?|+')). (3.31) 

3. We use compactness to prove that the sequence (g", u") converges, up to extraction, 
to a solution of (1.6). 

3.2.1 Step 1: Friedrichs approximation 

In order to construct approximated solutions of system (1.6) we shall use the classical 
Friedrichs approximation where we define the frequency truncation operator J„ by: 

for all n G N and for ah g G L^R''), J^g = T'^ (li<|.|<„(05(e)) , 



-^<I?I<"V 
and we define the following approximated system: 

dtqn + Jn{JnUn ' ^ JnQn) + J-rM'^ + Jngn)div J„U„ j = 0, 

dtUn + Jn{JnUn ' ^ JnUn) - AJnUn - 2 J„ f Z:'( J„ti„).V(ln(l + .InQn)) 

+ Jn{v{G{l + Jnqn))) =0, 
, {Qn,'Un)t=0 = {JnQO, JnUo). 

we recall the operator A is defined by: Au = An + Vdivu. 

We can easily check that it is an ordinary differential equation in L^^ x (L^J , where 
L'^ = {u (z L'^{W^), JnU = u]. Then for every n G N, by Cauchy-Lipschitz theorem 
there exists a unique maximal solution in the space C^([0,Tj*[, L^J and this system can 
be rewritten into: 



dtQn + Jn{Un ' ^Qn) + Jn (^(1 + gn)divn„j = 0, 

dtUn + Jn{un ' Vu„) - Aun - 2 J„ f i:'(n„).V(ln(l + qn))j (3.32) 

+ J„(v(G(l + g„))) =0. 

3.2.2 Step 2: Uniform estimates 

In the sequel, we will split n„ into the solution of a linear system with initial data JnUo, 
and the discrepancy to that solution. More precisely, we define by uj the solution of the 
following heat equation: 

dtul -Aul = 

{UL)/t=o = JnUo- 
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We now set Un = Un — u^- Obviously, the definition of Un leads to the following system: 



dtqn + Jniiun + ul).Vqn] + J„f (1 + qn)div {un + n^)) = 0, 

dtUn - AUn - Vdivu„ + JnUun + ul).VUn] + Jnlun-^ulj + Jn(ul.\7ul' 

- 2J„(D(n„ +O.Vln(l + (7„)) + VJ„(g(1 + g„)j = 0, 

, {qn,Un)t=0 = iJnqO,0)- 

(3.34) 
We would like to obtain uniform estimates on {qn,Un) in the space Et (see 3.31). Before 

doing this, let us recall that thanks to proposition 4 and as JnUo uniformly belongs (for 

• -—1 
all n) to -Bqo^i n B22 we obtain that for all T > 0: 

<C{\\uo\\.a_, . ^ +11/11 ._, . ). (3.35) 

In particular by Besov embedding, we can remark that Vu^ belongs in L^(L°°), this 
property will be crucial in the sequel in order to estimate u„- We would also point out 

that as div J„iio uniformly belongs (for all n) to i?! 1 we obtain that for all T > 0: 



||divu7|L dn +||divu?|| d,i <C(||uo|| d_, +,,„,, 

(3.36) 
Thanks to the apriori estimates from lemma 2 (as Jnqn = Qn and JnUn = Un the proof 
of this lemma, which is based on L^-scalar products, remains true) there exists r/ > 
and a time T > (all of them independant of n) such that for any n G N and any 
te [0,mm{T*,T)], 

\qnL .^ <e3^Ml + ||go||.4 ) -1, 

^tB2^i V ^2^1/ (3.37) 

|'Un|U_^ .d_i + ||u„|| ^ .d+i < 2r?. 

^?°-^2,l ^t-^2,1 

From this we deduce that the L^-norm of {qn,Un) is bounded. As J^ is the truncation 
operator in {^ G R'^, ^ < |^| < n} the bound blows up as n goes to infinity, but all that 
is important is that it implies (by contradiction) that for all n, the maximal lifespan 

n — • 

3.2.3 Step 3: Time derivatives 

Once the uniform time T is obtained the rest of the method is very classical. Using the 
previous uniform estimates to bound the time derivatives of the approximated solutions, 
we obtain that: 

Lemma 3 With the same notations, {dtqn)n is (uniformly in n) bounded in Li^Bi^i 
and {dtUn)n is bounded in L|,i?2^^ ^ + L'^B!^^ and then in Lj.{B2i + B21 
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d_3 
2 ' 



This result allows to get that: 

* Qn — 9n(0) is (uniformly in n) bounded in CtB!^^ Pi Cr^B^i , 

• Un is bounded in CxBi^i r\C:^{B2i +-^2^1 ^)- 

3.2.4 Step 4: compactness and convergence 

This part is also classical and we refer for example to [18] (chapter 10) for details: using 
the previous result and the Ascoli theorem, we can extract a subsequence that weakly 
converges towards some couple (<?,«), which is proved to be a solution of the original 
system and to satisfy the energy estimates. This concludes the existence part of the 
theorem. 

4 Uniqueness 

Once more, system {SW) is very close to {NSC) and the uniqueness is dealt the same 
way except, obviously, that here the external force terms that have to be plugged into 
the apriori estimates are different from the ones in {NSC) and we will focus on it in this 
section. As for {NCS) we will have (due to endpoints in the paradifferential remainder) 
to treat separately the cases d = 2 and d > 3. The second difficulty is that, as we present 
a local result for large data, we will have to make frequency cut-off (as in [2]) in order to 
bound some external force terms. 

Let us introduce for s G M the following space: 

ES) = Lr{Bl2' n B^\) n lI{bI^' n bI^,). 

Theorem 3 Let d>2 and assume that {qi,Ui) (i £ {l,2}j are two solutions of {SW) 
with the same initial data on the same interval [0, T] and both belonging to the space 
Ed{T). Then {qi,ui) = {q2,U2) on [0,T]. 

Proof: for i £ {1,2}, let us introduce Ui = Ui — ul (see (3.25) for the definition of 
ul), then {qi,Ui) satisfy the system: 

' dtqi + {ui + UL)NJqi + (1 + gj)div (uj + ul) = 0, 
dtUi - Aui - Vdivuj + {ui + UL).Vui + Ui.Vui + ulN/ul 

- 2D{u, + UL).v(ln(l + g,)) + v(g(1 + qi)) = 0, 

and if we denote by Sq = qi — q2 and 5u = ui — U2^ then {6q^Su) satisfy the following 
system: 

f dt5q + {ul + U2)-y5q = 5Fi + 5F2 + ^Fg, 

y dtSu - A5u - Vdiv 6u + {ul + ui).V5u + 5uNJ{ul + U2) = 5Gi + 5G2 + ^Gs, 

(4.38) 
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with: 



6Fi = —5u.Vqi, 

6F2 = — (1 + qi)div6u, 

SFs = -6q.div {u2 + ul), 

6G1 = 2D{uL + ui).V ( ln(l + gi) - ln(l + 92 ! 

6G2 = -Di6u).V\nil + q2), 

SGs = -V(G{l + qi)-G{l + q2)). 



4.1 The case d>3 

We wish to prove (as for (NSC)) the uniqueness in the following space: 



i _i / .a _9 . d \ d 



Ft = CtBI^ X [CtBI^ n L^5|i 

Due to endpoint estimates for the paradifferential remainder, the case d = 2 has to be 
treated in a different space. We refer to the following section. As for the classical {NSC) 
system, we prove that {5q, 5u) G Ft (the proof is left to the reader and the computations 
are the same as the ones done in the following). 

Let us begin with 5q. As qi and §2 have the same initial data, doing the same 
computations as for the transport estimates (see [2] theorem 3.14) leads to: 



< 



/ ( (||div(nL + tZ2)||L- + \\yiu2+UL)\\ . 



S2%nL- 



r 00 d"^ 
^T -02,1 



_-^ + \\5F\\ . d_-^ dr, 



(4.39) 



where 6F = SFi + 5F2 + 6F3 . These terms are estimated thanks to the paraproduct and 
remainder estimates recalled in section 2 (see (2.11)): 

-—1 
• Thanks to the Bernstein lemma we have B^l "^-t- B21 so that: 

|<J^ilLd_i < \\TsuVqi\\^d_j^ + \\T^q-^5u\\^d_-, + \\R{5u,Vqi)\\^d-i, 



^2,1 



^2,1 



B^' 

-°2,1 



<\\5u\\L^\\Vqi\\d_i + \\Vqi\\^-i \\6u\\ d + \\Su\\ d \\Vqi\\ d^^ 

-°2,1 -°2,1 -"2,2 ^2,2 

< 11(^^*1160 + ||5m|| d llgill d < \\Su\\ d llgill d . (4.40) 



• Similarly we get that 



SF2\\d_, < l + lki 



B^ 

-°2,1 



. d ■ 
-^2,1 



^i^sll .4-1 < ¥q\\ .d_^\\uL+U2\\ .4+1 



^2,1 



^2,1 



B?r<^Bi, 
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then we obtain: 






^+1 ., IML .^^, + (l + \\qi\\.^)\\5u\\.AdT. (4.41) 



< / \\UL +U2\\ .d_ 

Jo \ B? 

Finally, thanks to the Gronwall estimate: 

pglL d_, <e S2_^"s^,i / l + ||gi|| d )||(5«|| d dr. (4.42) 

Remark 7 Note that this estimate on 6q is vaUd for all d > 2. 

Concerning the velocity, using the a priori estimate for the transport-diffusion equa- 
tion provided in the present article, we can write that (we recall that ui and U2 have the 
same initial data): 

C/o (l|V«i+V«i.ll d +||V«2+V«l|| d )dT 

\\Sun .d_^ + \\5u\\ .d<e ^2.2^^^,! ^2,2^5^.1 



■ r oo d"? t 1 b"2" 

-^t -^2,1 ^t-^2,1 



X /" pGi + ^Ga + ^Gsll d_,,dT. (4.43) 

^0 ^2^1 



The last term is dealt the usual way: 



{SGsW .d_2 < C{\\qi\\Loo,\\q2\\L°°)U + \\qi\\ .d + \\q2\\ .d j\\Sq\\ .d_, 



<C(||go||..)||<^<?||.^-i. (4.44) 

BJ', B„^, 



Without surprise, the first term is estimated by: 



\\5Gi\\.d,<2\\D{uL+ui)\\d ||V ln(l + (7i)-ln(l + g2)n, .. 2, 

Bl^ Bl^r^&Ui ^ ' Bl^ 

<C(lko||.|)||^L+ni||.d+, \\bq\\.d (4.45) 

b}, B^l nfii 1 B,?, 

d o 

Note that here, after using the injection B^\^ ^-s- B21 , we needed d — 2 > in 
the remainder. In the case d = 2, this term will have to be dealt differently in the 
following subsection. 

As in (3.21), we have to decompose 6G2 into two parts: 

SG2 = -Di6u).V (ln{l + ^2) - ln(l + Smq2)) - D(<5n).Vln(l + ^^ga) = Ri + R2- 
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Here the second paraproduct for R2 requires that —1 + a < 0. The remainders 
require that d — 2 > 0). We obtain that: 



l-Rlll . d_2 < \\Su\\ .4 C'(lk2||L°°) 1+ 1^211 .4 Ik2 - Smq2\\ d 

^2,1 ^2,1 ^2,1 ^2,1 

< C'dlg'oll . d )||g2 - Smq2\\ .d \\5u\\ . d 
^2,1 ^2,1 ^2,1 



< 



C(lko||.4)(lko-5™go||.4 +(e^^(*)-l))||<5«||.4 , (4.46) 
and 



n2 ■ V ■02 / R- 



||i?2|| .4-2 < ll<^^^ll .4-cl|ln(l + 5'm^2)|| .4+„ < \\5u\\ d_^C{\\Smq2\\L°-)\\Smq2\\ . 4+„ 
^2^1 ^2^1 Bl^ Bl^ Bl, 

<C{\\qo\\d)2^n\m.d^- (4.47) 

R^ R^ 

The first term is smah if m is large enough and T smah enough, and the second 
term introduces a nonnegative power of t. 

Finally we have: 

C/o(l|V«i|| d +||V«2|| d +I|V«l|| d )dT 

M(-M M(-M ^ ^ R^ R^ R^ n R*^ / 

\\Sun d_^ + \\5u\\ d <e 2.1 2,1 S22ns^_j 

roo r7 ^ r 1 67 

-^t -°2,1 -^t-°2,l 

xC{\\qo\\.d) f[{l + \\Vui+VuL\\.d )||5(?||.d_,+2-™||<5ll||.d_ 

C/o*(||V«2|| d +||V«d| d )rfr \ 

+ (\\qo-Smqo\\ d +{e ""h ^h^^W - 1) ) ||,5n|| d dr. (4.48) 

^ Bli ' Bl^J 

Introducing cq a constant only depending on \\qo\\ . d and: 

-^2,1 

^(i) = / fllVuill d +IIVU2II d +||VttL|| d \dT. 

L V R^ "r^ r^ nRO / 

•/O ^2,1 -°2,1 ^2,2' '^00,1 



\\5uL .d_, + \\du\\ a <Coe^''^'U{t + Vm\6q\l . d_, 

-^t -'^2,1 ^t-'^2,l \ ^t ^2,1 

+ 2-™||<5n|| ^.d_ + (||(?o-5™go||.^ +(e^^W-l))||5u|| ^.d ) 

^tB2,l ^ ^2,1 ^ ^4-^2,1/ 



(4.49) 



If /3(t) = ||u||_ .d_i + ||ti|| . d^i, as in section 3.1 we have II (5tt 1 1 . d_^ < ta /3(t), then 

as, 

||5(?|U^.d_, <coe^^W||,5n|| ^.d , (4.50) 
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we obtain: 

m < coe^^^W (t + V{t) + 2"™tt + ||go - S„,qo\\ . . + (e^^W - l))/3(t). 

When f] G]0, 1] satisfies: 

2cor?e2^'' < \ (4.51) 

and m is chosen such that ||(7o — 'S'm'Zoll d < ri. Then if T is small enough so that: 

fn (||Vlir|| d +||divUr|| d ]dT<'n, 

° ^ ^^l^ni^So.i BlJ - (4.52) 

T + V{T) + (e^nT) _ 1) + 2"'"rt < r/. 

then we have, for all t € [0, T], 

< /3(t) < 2cor,e2^''/3(t) < ^. 
So /3(t) = for all t G [0, T] and the same goes for 5q, which proves the uniqueness on 

[o,r]. 

Remark 8 Note that these conditions are iniphed by the ones from the apriori estimates. 
To end the proof when T is not small, let us introduce (as in [18], section 10.2.4) the set: 

l'¥{t G [0,T]/iq,it'),Mt')) = iq2it'),Mt')), Vt' G [0,t]}. 

This is a nonempty closed subset of [0, T]. Using the same method as above allows to 
prove it is also open and then / = [0, T]. 

4.2 The case d = 2 

In this case, the estimates on 5q remain correct, but the paradifFerential remainders, when 
estimating the external forces terms in the velocity equation, are modified. Indeed, in the 
case d = 2 we reach the following endpoint where for all 1/pi + l/p2 = 1 = 1/ri + l/r2: 



Ds iiy II R^^ 



\\Rif,9)\\B0 <C\\f\\B^ \\9 

i,oo 

Estimate (4.43) is then replaced by 



|._ll , lU-ll ^ C'/(,M||V«i+V«i||^i ^^0 +||V«2+V«i,||^i ^^0 J'iT 

\Ou\\roor,-l +\\Ou\\rlr,l <e V 2,2 oo.l 2,2 oo,!/ 



i-t 
X 
'0 



/ \\5Gi + 6G2 + 5G3II 0-1 dT, (4.53) 
Jo ^'°° 



with 

ll'^G'sll^-i < C(||^i||Loo,||^2||L-)(l + |kl|lBi + h2\\Bi)\\Sq\\B0 

2,00 \ 2,1 2,1 / 2,00 

<C{\\qohiJ\\5q\\^o^, (4.54) 
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and (here we reach the endpoint d — 2 = in the remainder) 
\\SGi\\^-i = 2||Z)(ni +TIi).v(ln(l + q^) - ln(l + (^2)) h-^ 

2,00 \ / 2,CXD 

<2(||rz5V||^ 1 +||TvZ?||^-i +\\R{D,V)\\^i) 

2,0c 2,0c 1,00 

<2(||Z)||l.o||V||o-1 +||V||o-2 \\D\\r,l +\\D\\r,l ||V||o-l) 

<C(lkobi )||nL + ni||^2 nfii \M\bo- (4.55) 

2,1 2,2 oc,l 2,1 

Concerning the last term, with the same decomposition, 6G2 = Ri + R2 and when we 
choose some a g]0, 1[ (for R2) we obtain that: 



||^i||d-i < ||5li|Ui C(||g'2||L°°) 1 + 1^21101 Ik2 - 5'm, 92 II Di 
-°2,oo -°2,oo \ -°2,1/ -°2,1 

<C(|ko||«i )(lko-5„^go|lBi +(e^''^*^-l))ll'5^IUi , (4.56) 

J,l \ z.l / J,oc 

and 



p2b-i < ||(5n||^i-.||ln(l + 5„g2)bi+" < C-dltzobi )2"™||5tl||^i-.. (4.57) 

2,00 2,00 2,1 2,1 2,00 

As in the previous section, we coUect the estimates and obtain: 

||5^|IZc„so_^<coe^^(*)||<5n||^.^i_^, 

and with obvious notations, 



m < coe^^W ( j^ (l+V^'(T))||5g||j^^^.,^dr+(2°™tt + ||go-5„^goll4i^^ + (e^''(*)-l))/3(i) 
under the conditions from the previous section, we get: 

/3(t) <coe^^Wn\l + y'(r))||5n||^i^.^drj +i/5(t). 
Thanks to Proposition 3, we can use the fohowing logarithmic estimates {d = 2) 



\\5u\\ , <C\\5u\\ a log(e + 



\\6u\\ d , + ||5u 



iR^ "riR^ 



^t-^2,00 ^t-S; 



+1 , 



2.00 



t 2,00 



i}-B'?i l]b} V ||(5tt|| 

t 2,1 t 2,00 II r 1 j 

As 6u = ui — U2, we can write: 

||<5n||M.o +||<J^||.i.2 <W{t) = Wiit) + W2it), 

t z,oo t z,oo 

with Wi(t) = Iteillrino + Il^i'illrin2 . This function is bounded on [0, Tl and the 
estimates turn into: 

m < Ct J\l + F'(r))/3(T) log (e + ^^)dr. 
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As we have, 

"1 dr _ 

rlog(e + m)"°°' 

The Osgood lemma allows us to conclude that /3(t) = for all t G [0, T] (we refer for 
example to [2], section 3.1.1). Then the density fluctuation is also zero on this intervall. 
Then the conclusion is the same as in the case d > 3. 

5 Global well-posedness 

In this section we are interested in proving the global well-posedness of (1.1) when we 

assume smallness on the initial data. The proof follows the same lines than in the sections 

3 and 4. The only difficulty consists in getting damped effects on the density in order 

to deal with the pressure in the remainder terms. To do this we have just to use the 

estimates in Besov spaces from [14] or [6] on the following linear system associated to 

(1.1): 

dtq + V ■ Vq + divn = F, 

(5.58 
dtu + v-Vu- Au + Vq = G, 

There, they exhibit the parabolic smoothing effect on u and on the low frequencies of 
q, and a damping effect on the high frequencies of q. To do this, the authors need to 
introduce a paralinearisation in order to deal with the convection terms u ■ Vq. More 
precisely they obtain the following proposition: 

Proposition 5 Let {q,u) a solution of the system (5.58) on [0, r[ , 1 — | < s < 1 + ^ 
and V{t) = L ||Vf (t)|| d dr. We have then the foUowing estimate for any T > 0: 

\\{q,u)\\^^ ~s-i,s ^^ .^_^ +\\{q,u)\\^_^ ~s+i,s ^^ .^^^ 

<Ce''''^'\\\{qo,no)\\.s-us ^ + T e-^^W||(F, G)(t)||..-m . ^dr), 

~S1,S2 

where C depends only on N and s and B2 1 denotes the hybrid Besov space with 
regularity si for low frequencies and S2 for high frequencies (we refer to [14] or [6] for 
details). 

Using this on u, the rest of the proof simply consists in treating the remainder terms as 
in section 3, for the uniqueness the method follows the same approach as in section 4. 

6 Proof of theorem 2 

Our method from section 3 and 4 may be adapted to the study of incompressible density 
dependent Navier-Stokes equations. This is just a matter of replacing the parabolic model 
below by a nonstationary Stokes system. More precisely we define ul as the solution of 
the following system: 

dtUL - AuL + VUl = 0, 

divuL = 0, (6.59) 

(ul) t=o = uo- 
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In the same way than in section 3, we are searching solution of the form u = ul + u with: 

' dtq + {u + UL).Vq + (1 + g')div {u + ul) = 0, 
dtu — Au+ (u + ul).Vu + u.Vui + u^.Vul 

-2D{u + UL).v(ln{l + q)) + -^V{U + Ul) = 0, ^^'^^^ 

V y l + q 

div-u = 0. 

By applying the operator curl to the momentum equation, we obtain that: 
' dtq +{u + UL).Vq + (1 + g)div {u + ul) = 0, 



dtcurlu — Acurlu + curll u + ul)SJu + uSJul + ul.Vul 



2D{u + UL).V(ln{l + q)) )+V(^) : V(n + Hl) = 0, 



(6.61) 



, divu = 0. 



~ • #-1 
By following the same idea as in section 3, we are able to estimate u in L^{B2i ) n 

, . . —-L^ ■ ■ • — 9 - ■ ■ — 

L}p{B2i ) by proving estimate on curlu in L^{B2i ) n L}p{B2i), in order to deal with 
the pressure n it is sufiicent to adapt the idea of [1, 23] as IT verifies an elliptic equation. 

Aknow^Iedgements: The authors would like to thank Raphael Danchin and Pierre 
Germain for many fruitful discussions. 

References 

[1] H. Abidi and M. Paicu. Equation de Navier-Stokes avec densite et viscosite variables 
dans I'espace critique. Annales de I'institut Fourier, 57 no. 3 (2007), p. 883-917. 

[2] H. Bahouri, J.-Y. Chemin and R. Danchin. Fourier analysis and nonlinear partial 
differential equations, Grundlehren der mathematischen Wissenschaften, 343, Springer 
Verlag, 2011. 

[3] J.-M. Bony, Calcul symbolique et propagation des singularites pour les equations aux 
derivees partielles non lineaires, Annales Scientifiques de I'ecole Normale Superieure. 
14 (1981) 209-246. 

[4] J. Bourgain et N. Pavlovic. Ill-posedness of the Navier-Stokes equation in a critical 
space in 3D. J. Fund. Anal. 255 (2008), no. 9, 2233-2247. 

[5] D. Bresch and B. Desjardins, Sur un modele de Saint- Venant visqueux et sa limite 
quasi-geostrophique. C. R. Math. Acad. Sci. Paris, 335(12): 1079-1084, 2002. 

[6] F. Charve and R. Danchin, A global existence result for the compressible Navier- 
Stokes equations in the critical L^ framework, Archive for Rational Mechanics and 
Analysis, 198(1), 2010, 233-271. 



24 



[7] J.-Y. Chemin, Theoremes d'unicite pour le systeme de Navier-Stokes tridimensionnel, 
J.d'Analyse Math. 77 (1999) 27-50. 

[8] J.-Y. Chemin, About Navier-Stokes system, Prepublication du Laboratoire d' Analyse 
Numerique de Paris 6, R96023 (1996). 

[9] J.-Y. Chemin and N. Lerner, Flot de champs de vecteurs non hpschitziens et equations 
de Navier-Stokes, J. Differential Equations, 121 (1992) 314-328. 

[10] Q. Chen, C. Miao and Z. Zhang, Global well-posedness for the compressible Navier- 
Stokes equations with the highly oscillating initial velocity, Communications on Pure 
and Applied Mathematics,\o\ume 63(2010) 1173-1224. 

[11] Q. Chen, C. Miao and Z. Zhang, Well-posedness in critical spaces for the compress- 
ible Navier-Stokes equations with density dependent viscosities, arXiv :0811.4215vl 
November 2008. 

[12] R. Danchin, Fourier analysis method for PDE's, Preprint, Novembre 2005. 

[13] R. Danchin, Local Theory in critical Spaces for Compressible Viscous and Heat- 
Conductive Gases, Communication in Partial Differential Equations, 26 (78), 1183- 
1233, (2001). 

[14] R. Danchin, Global existence in critical spaces for compressible Navier-Stokes equa- 
tions, Inventiones Mathematicae, 141, pages 579-614 (2000). 

[15] R. Danchin, Global existence in critical spaces for compressible viscous and heat- 
conductive gases, Archiv for Rational Mechanics and Analysis, 160, pages 1-39 (2001). 

[16] R. Danchin. Zero Mach number limit in critical spaces for compressible Navier-Stokes 
equations, Annales Scientifiques de I'Ecole Normale Superieure, 35, pages 27-75 (2002). 

[17] R. Danchin, On the uniqueness in critical spaces for compressible Navier-Stokes 
equations. NoDEA Nonlinear Differentiel Equations Appl, 12(1):111-128, 2005. 

[18] R. Danchin, Well-Posedness in critical spaces for barotropic viscous fluids with truly 
not constant density. Communications in Partial Differential Equations, 32:9,1373- 
1397. 

[19] H. Fujita and T. Kato. On the Navier-Stokes initial value problem I, Archive for 
Rational Mechanics and Analysis, 16, (1964), 269-315. 

[20] B. Haspot, Cauchy problem for Navier-Stokes system with a specific term of capil- 
larity, M3AS, Vol. 20, No. 7 (2010) 1-39. 

[21] B. Haspot, Cauchy problem for viscous shallow water equations with a term of 
capillarity ,Hyperholic problems: theory, numerics and applications, 625-634, Proc. 
Sympos. Appl. Math., 67, Part 2, Amer. Math. Soc, Providence, RI, 2009. 

[22] B. Haspot, Existence of global strong solutions in critical spaces for barotropic vis- 
cous fluids, accepted in Archive for Rational Mechanics and Analysis 



25 



[23] B. Haspot, Local well-posedness for density-dependent incompressible fluids with 
non-Lipschitz velocity, accepted inAnnales de I'Institut Fourier. 

[24] B. Haspot, Well-posedness in critical spaces for the system of compressible Navier- 
Stokes in larger spaces, accepted in Journal of Differential Equations. 

[25] Akitaka Matsumura and Takaaki Nishida. The initial value problem for the equations 
of motion of compressible viscous and heat-conductive gases. J. Math. Kyoto Univ., 
20(1): 67-104, 1980. 

[26] Akitaka Matsumura and Takaaki Nishida. The initial value problem for the equations 
of motion of compressible viscous and heat-conductive fluids. Proa. .Japan Acad. Ser. 
A Math. Sci, 55(9):337-342, 1979. 

[27] A. Mellet and A. Vasseur, On the barotropic compressible Navier-Stokes equations. 
Comm. Partial Differential Equations 32 (2007), no. 1-3, 431-452. 

[28] Y. Meyer. Wavelets, paraproducts, and Navier-Stokes equation. In Current develop- 
ments in mathematics, 1996 (Cambridge, MA), page 105-212. Int. Press, Boston, MA, 
1997. 

[29] J. Nash. Le probleme de Cauchy pour les equations differentielles d'un fluide general. 
Bull. Soc. Math. France, 90: 487-497, 1962. 



26 



